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Abstract. We consider a relativistic hydrogenic atom in a strong magnetic field. The 
ground state level depends on the strength of the magnetic field and reaches the lower 
end of the spectral gap of the Dirac-Coulomb operator for a certain critical value, the 
critical magnetic Geld. We also define a critical magnetic field in a Landau level ansatz. 

In both cases, when the charge Z of the nucleus is not too small, these critical 
magnetic fields are huge when measured in Tesla, but not so big when the equation 
is written in dimensionless form. When computed in the Landau level ansatz, orders 
of magnitude of the critical field are correct, as well as the dependence in Z. The 
computed value is however significantly too big for a large Z , and the wave function is 
not well approximated. Hence, accurate numerical computations involving the Dirac 
equation cannot systematically rely on the Landau level ansatz. 

Our approach is based on a scaling property. The critical magnetic field is 
characterized in terms of an equivalent eigenvalue problem. This is our main analytical 
result, and also the starting point of our numerical scheme. 
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1. Introduction 



It is widely accepted in physics that eigenstates in strong magnetic fields can be well 
approximated in a Landau level ansatz. In this paper, we investigate the validity of such 
an ansatz in the case of the magnetic Dirac-Coulomb equation for large magnetic fields. 
More specifically, we study the critical threshold for the magnetic field defined as the 
smallest value of the magnetic field for which the lowest eigenvalue (the ground state) 
in the gap of the Dirac operator reaches its lower end. After a reformulation in terms 
of an equivalent minimization problem, we compare the values of the critical magnetic 
field with and without ansatz, both from an analytical and a numerical point of view. 

Our main theorem is stated in Section 2. It characterizes the critical magnetic 
field as a function of the lowest energy for an explicit eigenvalue problem. Section 3 
will be devoted to its proof. A Landau level ansatz is then defined in Section 4 and 
some comparison results for the critical magnetic field, with and without ansatz, are 
given there. A numerical method based on Theorem 1 has been implemented in both 
cases. Results are given in Section 5 and compared with earlier works, see (Schliiter 
et al. 1985). They show that for large coupling constants (in the electrostatic field), the 
critical threshold is well below the critical threshold in the Landau ansatz. 

Critical magnetic fields arc huge and can eventually be encountered only in some 
extreme situations like magnctars, which arc neutron stars with intense magnetic fields, 
see (Kouveliotou et al. 1999). This is the only known domain of physics for which our 
computations might eventually be relevant, see (Dolbeault et al. 2007) for a discussion. 
In the Dirac-Coulomb model, the value of the critical magnetic field only provides an 
order of magnitude of the field strength for which pair creation could eventually occur. 
Such a phenomenon should of course be studied in a full QED framework and the 
computations done using the Dirac-Coulomb operator are only an indication on the 
scales that should be taken into account. See (Pickl 2005) for more details. For a review 
on large magnetic fields in physics, see (Duncan 2000). 



2. Main results 

The magnetic Dirac operator with Coulomb potential iy/\x\ can be written as : 

, l-iy/\x\ -i(7-(V-zA) 
' -ia-{V-iA) -I-u/\x\ ' 

where A is a magnetic potential corresponding to B, and I and cxfc are respectively the 
identity and the Pauli matrices 
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Let B = (0, 0, B) be a constant magnetic field and the associated magnetic potential. 
For any x — {xi,X2,X3) &M.^, define 



Pb -.^ -ia ■ {W -iABix)) , Ab(x) := — 



V 



Xi 





(3) 



and consider the functional 



|Pr</)P o 



J[cl>,X,u,B]:= I |-^l|^ + (l-A)H^-^|0n d^x (4) 

on the set of admissible functions A{i', B) := {0 e C^(M^) : ||0||^2(r3) = 1 , A i-^ 
J[(f), A, B] changes sign in (—1, +oo)}. The essential spectrum of Hb is M \ (—1, 1). 
By (Dolbeault et al. 2007, Theorem 1), the smallest eigenvalue in (—1, 1) of Hb is 

Xi{u,B):= inf \[<P,u,B] (5) 

where A = A[(6, u. B] is cither the unique solution to J[(f), A, z/, 5] = if G Alu, B), or 
\[(f), u, B] = —1 if J[(f), — l,u,B] < 0. Also see (Dolbeault et al. 2007) for the relation of 
the two-components spinors in ^(t^, B) with the four-components spinors and the action 
of Hb on them. 

The critical magnetic field is defined by 



B{u) := inf |S > : limmf Ai(i/, b) ^ -ij . (6) 

Define the auxiliary functional 

Sb,M-- [ ^-^iPscPl'd'x- [ ^\cf>\'d'x, (7) 
Jr3 V \x\ 

that is SbA<P\+'^ UWl^m?) = J[<P. -1, B] . The scaling 0^ := E^/^ (5^2 x) preserves 
the LP' norm, and yields 

SbAM = ^^iA<P\ ■ (8) 

We define 

Recall that \i{v^B) is characterized as an eigenvalue of Hb only as long as it takes 
values in (—1, 1). If we could take the limit B B{u), we would formally get that 
— 1 = Ai(z/, i?(z/)), which, still formally, amounts to info^(/,ec^(R3) J[0, — 1, i/, 5] = 0. It 
is therefore natural to expect that ^J Biy) fx{v) +2 = 0. Proving this is the purpose of 
our main result. 



Theorem 1 For all v e (0, 1), ^i{y) is negative, finite 

4 

and B{u) is a continuous, monotone decreasing function of v on (0, 1). 



= tt::^ (10) 
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3.1. Preliminary results 



Lemma 2 On the interval (0, 1), the function v 
decreasing and takes only negative real values. 



is continuous, monotone 



Proof. The monotonicity of the function //(i/) is a consequence of its definition. On the 
other hand, the functional 

,,2 



(i/, 0) 1-^ PSl^iy[ 



\x\ |Pi(/)Pd^X - 



m3 FI 



\(l)\'d-'x 



is a concave, bounded function oi u e (0, 1), for any e C^(R^), and so is its infimum 
with respect to 0. A bounded concave function is continuous. 
Next, consider the function 



271 



\ J 



\J X 



a;i,X2,X3) e (12) 
1. Note that 



for some / G C^(M, M) such that / = 1 for |a;| < 5, 5 > 0, and ..^ ii^2(x+) 
(p G Ker(PB + ^ (73 d^-^) and so, Pb0 = — i cts (9i:30. Moreover, the function is normahzed 
in L^(R^). Same computations as in (Dolbeault et ah 2007, proof of Proposition 6) show 
that 



SBM<— + C2V-C^v\ogB, 



(13) 



where Cj, i = 1, 2, 3, are positive constants which depend only on /. For B > 1 large 
enough, £bA<^] < 0- With </)i/b(.x) = B-^l'' (t> {B-^^ x) , E^,,[<P^,b\ = B-^" EbM < 
by (8). This proves that ^{v) is negative. 



Lemma 3 For any a, 6 G R*^ and any £ > 0, 



\a?> 



■ ■ l+£ e 
Proof. A simple computation shows that 



eil + e) 



|a + &P 1^ 

1 + £ £ 



(£0-6)2 > . 



□ 

(14) 

(15) 
□ 



3.2. Proof of Theorem 1 

Let us first prove that B{u) < Bi^u) := A/ By definition of ji{u) and (8), 

B{u) = sup |s > : inf {£b,M + 2 1101112(^3)) > o| . (16) 

By definition of B(u), B < B(u) means that \i{u,B) > —1. Since Ai(z/, 5) is the 
infimum of X[(l), u, B] with respect to G C^(R^) and J[0, A, u, B] is decreasing in A, 

£bA4>] + 2 ||0||i2(M3) > J[0, Ai(z/, B), i/, 5] > J[0, A[0, i/, 5], i/, = (17) 
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for all 0, so that B < B{v). This proves that B{v) < B{v). 

To prove the opposite inequality, let B — B^u) and consider a sequence (^'n)neN 
such that G (0,1^), lim^^oo^n — ^-i '■— ^i{^n,B) > —1 and lim„_»ooA" = —1. 
Let (f)n be the optimal function associated to A": J[d>n, A", z/„, B] = 0. 

We define a sequence of truncation functions (Xn)neN as follows. Consider first a 
nonnegative smooth radial function x on M"*" such that x = 1 on [0, 1], < x < 1 and 
X = on [2, oo). Then we set Xn^x) := x(|,t|//?„) for some increasing sequence {Rnjnevi 
such that lim„_>oo Rn = oo. By applying Lemma 3 to 

Pb <t>n = {.Pb 4>n) Xn + l ' (^0 Xn) <Pn\ , (18) 
=a =b 

where (pn-=4>nXn, wc get 

|p^0„r>KMO2^_K,^MM, (19) 

f ^~ 

for some sequence (£n)neN of positive numbers, to be fixed. Hence, using the fact that 
< Xn < 1> we get 



/ 



l + A" + g 



1 + ^3 1 + A" + 1^ Jm3 1 + A" + 1^ 

The function 0„ is supported in the ball 5(0, 2 i?„): with /i„ := (l+£„)[2(l+A'*)i?„ + z/„], 

'^'^-1' £x>—l \x\\PB4>n?d''x. (21) 



1 + e„ Jkb 1 + A" + - /in 

We choose Sn and i?„ such that 

hm £n = , hm = oo and hm (1 + A") i?n = , (22) 

n— *oo n— »oo n-+oo 

SO that 

lim /In — ■ (23) 

n— *oo 

The function (PoXn) is supported in S(0, 2i?„) \ S(0, i?„): there exists constant k 
depending on ||x'||l°°(i,2) such that jPoXnP < i^Pn^ ^ ^ consequence. 



1 f |(Po Xn)0nP ^: 
£n 7183 1 + A" + ^ 



^a; < ^ ^- / d^x . (24) 

1^1 - £„i?„[(l + A-)i?„ + f]7M3' ' ^ ^ 

Moreover, 

^n/ ^-^d'x-uj \^d'x<^[ \4>n\'d'x. (25) 
Jm3 pI JkS \x\ Kn Jk3 

Thus, with rin ■— n/ {sn Rn ((1 + A") Rn + ^)) + Un/Rn, WC Can write 



= J[(f)n, A", //n, B]>—[ \x\ \Pb ^nT C^^^; - Vn I 

A*n Jm3 Jm 



M3 pI 



+ (l-A"-77„) / l^npd^x. (26) 

JIR3 
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Assume further that 

lim En Rn — oo , (27) 

SO that 1 — A" > r/n, — > as n — > oo. Using again the fact that < < 1, we get 



= J[0„, A'*, Un, B]>— \x\ \Pb (t>n? d^x - V, 



, [ M 
" / II 

Jm? fI 



+ (1 - A^ 



|2 j3 



d'^X . 



(28) 



Let Vn = Vn- We have obtained 



^/ \x\\PBKfd''x-v^\ \^d^x + 2l \4>n\'d^x 

Jr3 Jr3 \x\ 



< 



2- 




;i - A" - Vn) 



\(f)n\'^d^X 



(29) 



The left hand side is a decreasing function of i>„. Since hm„^ooi/„ = i/, for any i/' > i/, 
for n large enough, 



1 /' \x\\PB4^n\'d'x-:.' [ ^ " 



d^x + 2 I \k}'^d^x 



< 



\(pn\ d'^X . 



(30) 



We observe that by construction ^„ is non trivial. By homogeneity, one can even assume 



that 



= 1. Since lim„^oo J 'r (1 - A'^ - rj^) = 2, 



B{u') <B = B{u) ^ v' >v . (31) 
By Lemma 2, v' ^ B{i'') is continuous. This proves that B{u) < B{y). □ 

4. A Landau level ansatz 

In analogy with what is done in nonrelativistic quantum mechanics, we denote by first 
Landau level for a constant magnetic field of strength B, see (Dolbeault et al. 2007), 
the space of all functions which are linear combinations of the functions 

B 



[X2 



v^272^ 

where the coefficients depend only on X3, i.e., 

^(^) = /K^s) (t)e{xi, X2) . 



+ ^a;i)^e-^^'/^Q, £eN, s^ = xl + x 



2 

2 ' 



(32) 



(33) 



In this section, we shall restrict the functional £b,u to the first Landau level. In this 
framework, that we shall call the Landau level ansatz, we also define a critical field by 

(34) 



Bc{u) := inf <{ S > : liminf A^'(z/, h) = -I 
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Af (i/, B) := inf A[(/), u, B] . 

<i>£A{v,B) , n^(P=0 



(35) 



Here 11 is the projection of onto the first Landau level, and 11^ := I — IT. 

One can prove in the Landau level ansatz a result which is the exact counterpart 
of Theorem 1. For any v & (0, 1), if 

(36) 



then 



:= , inf ^i,z.[0] , 



The goal of this section is to compare //^(i^) with /x(i/) given by (9). By definition of 
these quantities, we have 

K^) < I^M . (38) 
With the notation s = ^/xl + and ^ = X3, if is in the first Landau level, then 



2 poo poo 

^iA(f>\ = y2~ hif^ dz-v I Qi fj dz , 
^ ^ Jo Jo 



where 



and 



at{z) := ( (j)e,-(l)e 



+00 ^2£+l g-sV2 



ds 



(39) 



(40) 



bAz) :-- 



(41) 



/•+00 

^'^r0,)^,(«, ,,,) = ^ s^^+' e-^'/Vs2 + z'' ds . 

A simple integration by parts shows that (0^, -F(r) 0^)^2(]g2 is increasing (resp. 
decreasing) in i whenever F{r) is increasing (resp. decreasing). Since and be only 
depend on \z\, iic{'^) is also achieved by functions which only depend on \z\ as well. It 



follows that 



- be jf dz-v ae fj dzj 

.if.o(i:/f)--f«.(i: 



f! dz 



Using the inequality 

d 



dz 



V^efe 



dz 



fi 



we find that 



1 POO POO 

SiM>- bo\ffdz-v aofdz, 
^ Jo Jo 



(42) 



(43) 



(44) 
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with / := \/^^yf- In other words, for our minimization purpose, it is sufficient to 
consider functions of the form 

H^) = m ^ ■ (45) 



/2¥ VO 

We observe that 



^ / \<l>\'d'x= r fdz (46) 



and 



- S,,M = - / hfdz-v afdz:^ , (47) 

^ ^ Jo Jo 

with a — ao, b — bo, i.e., 



b{z)^ Vs^ + z^se-' ds and a{z) ^ . ds . (48) 

Jo Jo V + 

The minimization problem in the Landau level ansatz is now reduced to 

/^/:M = inf^^#L. (49) 



It is a non trivial problem to estimate how close /i(z/) and f^c{^) are. Let 



/ -fdz-v a fdz (50) 

I' Jn a Jn 



and 



C[/]:=-/ hfdz-v -fdz, (51) 
^ Jo Jo 

with corresponding infima fJ'^ii') and //^(i^). 
Lemma 4 For an?/ e (0, 1), 

/i^(^^) </X£(i^) </x^(i^) . (52) 
Proof. This follows from 

6(z) > 4t (53) 

which in turn follows from Jensen's inequality, noting that s e"'*^/^ ds/ 2 tt is a probability 
measure. □ 



In (Dolbeault et al. 2007) it was proved that 
log|//Z(^)l ~ 

as z/ — > 0+. The methods in (Dolbeault et al. 2007) can be adapted to show that 

log|/^£(z^)l ~ 
as well, thus proving the following result. 

TT 

Lemma 5 With the above notations, lim u log|/X£(i/)| = 

i/-^o+ 2 



log\n^{u)\^ (54) 



log|/.+ (i.)|«-^ (55) 
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As 1/ — > 0+, this provides us with the leading asymptotics of the critical magnetic fields. 

log Biu) 

Theorem 6 With the above notations, lim = 1. 

i^^o+ log Bc{y) 

Remcirk. This result does not prove that ■§^^ converges to some finite limit as — > 0+ . 

Proof. By (38), we already know that ix{u) < fj.c{^)- To prove our result, we need 
an estimate of //(i/) from below. Since for i/ small, is very large according to 

(Dolbeault et al. 2007, Corollary 11), we can assume that B{u) > 1 for any u e (0, P) 
for some P > 0, so that. Hue (0, P), then Xi{u, 1) > — 1 and therefore, for all 0, 

SrA<P]>^M--= [ X r^l'^^f 1 . d'x- [ ^A<P\'d'x. (56) 

7r3 Al(i/, ij + i + 1^ \X\ 

For completeness, we sketch the main steps of the proof, which is similar to the one 
given in (Dolbeault et al. 2007). Since 



IS concave m x. 



1 + jp-,[</.] = sup — ^nr^ — • (^s) 



X l|V^llL2(1.3) + ||X||L2 

Obviously, 

G^it) ^ g-g) ,,,, 

sup -7— ^^T^ — > sup — — ^- — . (59) 



X + ||x||l2(r3) ni^=o ||0||l2(r3) + ||x||l2 

This, by (Dolbeault et al. 2007, Proposition 14), is bounded below by 

T l|n0|U2(K3) + ||n^0|U2(K3) + ||nxlU2(K3) 

Moreover, by (Dolbeault et al. 2007, Proposition 15), for u small enough, this is bounded 
below by 

||n0|U2(M3) + ||n^0|U2(M3) + ||nx|U2(M3) ' ^ ' 

where d{-) is a continuous function such that d{0) — -\/2- The inequality that leads 
to (61) displays the fact that being perpendicular to the lowest Landau level raises the 
energy. 

Again by concavity, for all 0, there is a unique x realizing 



T l|n0|U2(^3) + ||nxlU2(M3)-^w.3/2i 



, II =: A- ,3/. [0] . (62) 



Since A^^^g/j [0] < 1 < \/2, we finally have 

1 + ^40] >A:^,3/.[0] • (63) 
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(64) 



and therefore, for small enough, we have proved that 

l^-ci^ + 1^'^') < < fJcii^) , (65) 

where the upper estimate is given by Inequality (38). Since u — > 0+, u^^'^ becomes 
insignificant compared to v. The conclusion then holds by Theorem 1. □ 

5. Numerical results 

5.1. Computations in the Landau level ansatz 

To compute we minimize 'Ci^[/]/||/||^2('K+-) on the set of the solutions fx of 

(A + i/a(^))/ = 0, /(O) = 1 , /'(O) = . (66) 



b{z) ' b{z 
We notice that b'{z) = za{z), and, for any ^ > 0, 



+ z 



(67) 



Numerically, we use a shooting method and minimize (7(A,2;jnax) '■= |/A(-2max)P + 
|/I(^max)P for some ^^ax large enough. As ^r^ax ^ oo, the first minimum /i£(i/, ^^ax) of 
A H- >■ g{\, -2max) couvcrgcs to and thus determines A = fJ^ci^)- See Figure 1. 

m^c^ ~ 4.414 ■ 10^ be the numerical factor to obtain the critical field in 



Let b 



e h 



Tesla. Corresponding values are given in log^g scale. The minimum yU£(z/, 2;max) is found 
by dichotomy. Results computed with Mathematica are given in Table 1. 



v 


Z 


fJ-c 


Bc{y) 


logio(6B£W) 


0.409 


56. 


-0.0461591 


1877.35 


12.9184 


0.5 


68.52 


-0.0887408 


507.941 


12.3506 


0.598 


82. 


-0.14525 


189.596 


11.9227 


0.671 


92. 


-0.192837 


107.567 


11.6765 


0.9 


123.33 


-0.363773 


30.2274 


11.1252 


1 


137.037 


-0.445997 


20.1093 


10.9482 



Table 1. Numerical values found in the Landau level ansatz. 



There is no significant difference with the results that were found in (Schliiter 
et al. 1985). See Fig. 2. 
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Figure 1. Plot of log[l + ^(A, Zmax)] with, z^a.^ = 100, for v = 0.9. 
5.2. Computations in the unconstrained case 

Wc numerically compute B{v) in the general case, without ansatz. For this purpose, we 
discretize the minimization problem defining n{y) using B-spline functions of degree 1. 
They are defined on a logarithmic, variable step-size grid, in order to capture the 
behavior of the eigenfunctions near the singularity. We also use cylindrical symmetry to 
lower the dimension from 3 to 2. These two choices provide us with very sparse matrices, 
even if large. We use Matlab routines to calculate all integrals and the eigenvalues of 
the corresponding discretized matrices. The size of the computing domain is adapted 
as V varies. Results of these computations are shown in Table 2. 
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Figure 2. Left: values of the critical magnetic field in Tesla (logiQ scale) . Right: values 
in dimensionless units. Ground state levels in the Landau level ansatz correspond to 
the upper curve, while the ones obtained without approximation are given by the lower 
curve. Dots correspond to the values computed by Schliiter, Wietschorke Sz Greiner 
(Schliiter et al. 1985) in the Landau level ansatz. 



V 


Z 


Al 


B{v) 


logio(6i3(j.)) 


0.50 


68.5185 


-0.0874214 


523.389 


12.3637 


0.55 


75.3704 


-0.119458 


280.305 


12.0925 


0.60 


82.2222 


-0.153882 


168.922 


11.8725 


0.65 


89.0741 


-0.191037 


109.604 


11.6847 


0.70 


95.9259 


-0.231198 


74.833 


11.5189 


0.75 


102.778 


-0.274665 


53.0216 


11.3693 


0.80 


109.63 


-0.321875 


38.6087 


11.2315 


0.85 


116.481 


-0.373535 


28.668 


11.1022 


0.90 


123.333 


-0.430854 


21.5476 


10.9782 


0.95 


130.185 


-0.496005 


16.2588 


10.8559 


1.00 


137.037 


-0.573221 


12.1735 


10.7302 



Table 2. Results of the minimization method without symmetry ansatz. 
5.3. Discussion 

When dealing with the physics of magnetars, one is interested only in the order of 
magnitude of the critical magnetic field. With this goal in mind, the values computed 
in the Landau level ansatz are quite satisfactory. The corresponding values are given in 
the right column of Table 3, logiQ{b B{u)), where b is approximatively 4.414 • 10^ Tesla. 
See Fig. 2 (left). 

Except maybe in the limit — >■ 0, it is however clear from Theorem 1 that the 
equivalent eigenvalue problem (9) has nothing to do with its counterpart ^ciy) in the 
Landau level ansatz. What our computations show is that the values of the computed 
critical fields significantly differ, see Fig. 3, and that the shapes of the corresponding 
ground state do not have much in common, see Fig. 4. 
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4 ■ 




Figure 3. Ratio of the ground state levels computed in the Landau level ansatz versus 
ground state levels obtained by a computation without any symmetry ansatz. 




Figure 4. Plot of logio(10"^ + \())\'^) for v = 0.5 (top) and u = 0.9 (bottom). Left: 
Landau level ansatz. Right: computation without any symmetry constraint. Level 
lines are equidistributed. 



6. Conclusion 

Orders of magnitude of the critical magnetic field given by the Landau level ansatz are 
similar to the ones obtained without constraint. Qualitatively, the curve of the critical 
magnetic field in terms of the charge, is also well reproduced by the Landau level ansatz. 

However, the difference of the values of the critical magnetic fields is not small when 
comparison is done for values of u approaching 1. It turns out that values given by the 



Critical magnetic field in relativistic atomic physics 



14 



Landau level ansatz are 50 % higher than the ones found by a minimization approach 
without symmetry ansatz. Shapes of the density distributions also differ significantly. 
They are much more peaked around the singularity when computed unconstrained than 
in the Landau level ansatz. 

The Landau level ansatz, which is commonly accepted in non relativistic quantum 
mechanics as a good approximation for large magnetic fields, is a quite crude 
approximation for the computation of the critical magnetic field (that is the strength of 
the field at which the lowest eigenvalue in the gap reaches its lower end) in the Dirac- 
Coulomb model. Even for small values of which were out of reach in our numerical 
study, it is not clear that the Landau level ansatz gives the correct approximation at 
first order in terms of v. Hence, accurate numerical computations involving the Dirac 
equation cannot simply rely on the Landau level ansatz. 
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